Diffusive instability of a vortex in a rotating Bose gas 
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The dissipative dynamics of a vortex in a finite temperature trapped Bose-Einstein condensate are 
shown to be governed by a diffusive instability. In the weakly interacting regime we find a cross-over 
from instability to metastability of the vortex, determined only by the relative rotation rate of the 
vortex and thermal cloud. Even in the thermodynamically stable regime there is a finite lifetime 
arising from diffusive instability. We find the mean exit time for the vortex from the condensate, 
and show that near the critical frequency for thermodynamic stability the lifetime has a universal 
form which is independent of the condensate size and the vortex precession frequency. In the 
thermodynamically stable regime the steady-state two-time position correlation predicts 'thermal 
bunching', giving an increased probability of positional coincidence for short times. 

PACS numbers: 03.75.Kk, 67.40.-w, 32.80.Pj 



I. INTRODUCTION 

Since the achievement of atomic Bose-Einstein condensation (BEC) pj, the study of quantised superfluid vortices 
has attracted intense interest. The observation of vortices in BEC gave direct evidence of superfluidity [3- C3l - and led 
to studies of the dynamics of formation and decay of vortex lattices. Vortices also form a challenging testing ground 
of BEC theory which is renowned for its tractability, giving unique insights into quantum many-body physics. 

It is widely appreciated that there is an energy barrier to the penetration of a vortex into the interior of a BEC 3 . 
The experimental challenge of creating condensates containing large numbers of vortices was met by generating a 
thermal non-condensed fraction of the gas with significant angular momentum, which acts as a rotating environment for 
the condensate. Stimulated collisions between condensed and thermal atoms impart angular momentum to the BEC, 
leading to the formation of quantised vortices. Experimental methods include stirring with an elliptical perturbing 
potential p| , with a detuned laser [|| , and, alternatively, surmounting the barrier a priori by evaporatively cooling a 
rotating thermal cloud to quantum degeneracy . 

Vortex motion is strongly dependent on certain types of instability Q, and the role of instabilities in overcoming 
the energy barrier has been investigated theoretically and experimentally. There are two well-understood forms of 
instability in trapped Bose-Einstein condensates, dynamic and thermodynamic 0]. Dynamic instability arises when 
small perturbations of the condensate may grow with time, and is associated with an imaginary eigenfrequency of the 
linearised Bogoliubov-de Gennes (BdG) equations. Thermodynamic instability occurs when the system energy may 
reduce in the presence of dissipation, and is associated with negative eigenfrequencies. The aim of this paper is to 
show that the thermal motion of a vortex is strongly dependent on a new diffusive instability, and to elucidate the 
effect of thermal atoms on the vortex in a trapped rotating BEC. 

In this work we generalise the stochastic variational method developed by Duine et al. 0, to find an effective 
Lagrangian for the motion of a straight line vortex in a rotating Bose gas. We verify the validity of their approach for 
the regime treated in Ref . ^(j , and investigate the new physical processes revealed by the more general theory. 
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FIG. 1: Schematic of the system consisting of a partially condensed Bose gas in an oblate harmonic trap. The condensate 
contains a vortex which precesses at angular frequency uj p , and is immersed in a thermal cloud rotating at angular frequency 
f2. The circulation of the superfluid near the vortex core is indicated by v s . All quantities are defined with respect to the 
laboratory reference frame. 



II. THE SYSTEM AND VORTEX ACTION 
A. Approximate high temperature description 

There have been several theories developed to descibc the dynamics of a partially condensed Bo se g as. While 
exact theoretical methods, which in principle include all quantum coherences, are available [H [H Eg, practical 
calculations for nonlinear multimode systems using these methods are often severely limited by technical problems. 
The two approaches that interest us here are those of Stoof [3 E1j an d Gardiner 0, who have developed 
closely related classical field theories of Bose-Einstein condensation in the high-temperature regime. For a mode with 
characteristic frequency ui, the system is in the high temperature regime when the temperature satisfies 

Hlu < k B T < k B T c , (1) 

where T c is the temperature for Bose-Einstein condensation. The classical field approximation for the condensate 
mode is valid when it is macroscopically occupied, or when 1 <C N c . For our purposes the two theories are equivalent, 
so we will use the notation of Stoof to assist comparison with the paper of Duine et al |lfj |. The reader is referred 
to Ref. E3 f° r the details of the derivation of the SGPE for the scenario where the thermal cloud is at rest in the 
laboratory frame; the new physical feature we include here is the possible rotation of the thermal cloud. At the level 
of approximation we are using the rotation may be included by minor technical modifications of the theory, which 
nonetheless have important physical consequences for the motion of the vortex. 

The system of interest consists of a partially condensed Bose gas of atoms with mass m held in an axisymmetric 
oblate trap described by the potential 

^(x) = ^V+ y 2 )+^z 2 }, (2) 

where u> <C uj z . In the weak interaction regime the collisions are described by the two body T-matrix, T 2B — ATrh 2 a/m, 
which is determined by the S-wave scattering length a E3- The condensate is assumed to contain a single vortex 
which has a single quantum of circulation, and the thermal cloud at chemical potential /i is in a state of uniform 
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rotation about the z-axis at angular frequency Q with respect to the laboratory frame. This is shown schematically 
in Fig. H| 

For the general case where the thermal cloud is rotating, the theory of coupling between thermal cloud and conden- 
sate is most easily formulated in the frame of reference where the cloud is stationary. This leads to some modifications 
of the SGPE [l9j |. arising from the transformation into the rotating frame, in the form of terms which depend on the 
operator for the z-axis component of angular momentum, L z . In the frame of the thermal cloud, the essential change 
to the theory is to alter the single particle Hamiltonian for the Bosc gas to H sp — > H sp — QL Z [2J| • 

Incorporating the physical effect of the rotating cloud into the SGPE formalism, we take our starting point as the 
stochastic Gross-Pitaevskii equation 



dt \ 4 J { 2m 

where the fluctuations are completely determined by the non-zero correlator of the complex Gaussian noise amplitude 

(„•(*, tMx'.f)) = * n2 ^ (x '* ) g(x-x')g(t-0, (4) 

where the Keldysh self-energy arises from Bose-enhanced collisions of condensate atoms with thermal cloud atoms, 
and is given by @, OH E2 

RE*(x,t) = — 47n(T 2B ) 2 f [ f -^(2^^ - k 2 - k 3 ) 



,2tt) 3 J (2tt) 3 J (2tt) 3 
xS (ei - e 2 - e 3 )[/i(l + /a)(l + fa) - (1 + /O/a/a], (5) 

where the thermal cloud is described by the rotating frame Bose-Einstein distribution within the Hartree-Fock ap- 
proximation given by 

fi = [exp(-f3(ei-»N))-l]-\ (6) 
h 2 k? 



2m 



y cxt (x)+ 8 fin ■ (x x ki) + 2T iB |(0(x,i))| 2 . (7) 



form £1L£2,£ 



In this work we will treat the thermal cloud as time-invariant and spatially uniform. As |f2| increases the thermal 
cloud is c entrifug ally broadened, improving the validity of this description. We may then use the simpler approximate 

sK = -i48ma 2 (k B T) 2 

irh 3 { 1 

Equations |21 El and 0] describe a condensate driven by thermal fluctuations arising from Bose-enhanced collisions 
between condensate atoms and atoms in the rotating thermal cloud. The condensate is described by a modified 
Gross-Pitaevskii equation with appropriate damping and noise terms determined by the motion of the thermal cloud. 
It follows from the form of the noise and the damping terms that the SGPE automatically satisfies the fluctuation- 
dissipation theorem. We remark here that the interaction of the thermal noncondensate with the condensate leads to 
another process, in addition to the growth terms we are concerned with here, involving collisions which do not change 
the condensate occupation, the so-called scattering terms |l7j . However, as shown by Anglin and Zurek p3 |. in a first 
approximation we may neglect these terms since their net effect is to induce a small change in the effective damping 
and fluctuation terms in Eq. |3| 



B. Effective vortex action 



It has been shown in Ref. that the Wigner probability distribution for the condensate wavefunction evolving 
according to the Langevin field equation (jSJ) can be written as the functional integral 

W[4, 4*] = / c#M<f] exp -S eS [0, </>*] , (9) 



r, 
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where the effective action is 



dt' 



2m 



iJi- 



dt 1 



i + ^£ fc 



i/ ext + fii 2 + ^-r^|0(x,i 



(x,t) 



(10) 



In order to find the dynamics of a single vortex at finite temperature we introduce a Gaussian variational ansatz 
for the condensate wavefunction 



CM) 



^^l 2 q 2 q z [q 2 + u 2 x (t) + u 2 y (t)] 



{[x - u x {t)] + i[y - u y (t))) 



x exp 



x 2 + y 2 
2q 2 



2ql 



(11) 



which describes a single line vortex in an otherwise cylindrically symmetric Bose-Einstein condensate. In the nonin- 
teracting limit the ansatz reduces to the exact wavefunction for the condensate containing a vortex, since the length 
parameters reduce to the harmonic oscillator lengths determined by the trap. The time dependent variational param- 
eters u x (t),u y (t) are the Cartesian coordinates of the vortex, and the time independent variational parameters are 
the condensate size parameters in the radial (q) and axial (q z ) directions. Since the condensate can exchange particles 
with the thermal cloud it is crucial to also allow for the time dependence of the condensate occupation, N c (t) and 
phase (t). Substituting this ansatz into Eq. E]|we find the effective action 



S cn [N c ,9 ,u] 



dt' 



to 



N c (t') 



h9o(t') + Mc(*') - /i + -j (u x (t')uy(t') - u y (t')u x (t')) 



h 2 N c {t') 
q 2 

h 2 N c {t') 



N c {t') + ^^{^)-^)N c {t') 



'\l2 



[u x {t') + ui p (q, N c )u y (t') - iu x (t')] 



u v (t') - u p (q, N c )u x {t') - ju y (t')} 2 + 0{{ Ul /qf 



(12) 



where the damping rate, chemical potential, vortex precession frequency, and effective potential are respectively 

/3zS K 
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-h(u p (q,N c )-ty, 



(13) 



N c V(q, q z ,N c ) - h{w p (q, N c ) - O) xK ' yK ' 



dN, 



(14) 



v p (q,N c ) 



2mq 2 



2 ? 

muj q 

2h 



aN r h 



2nmq 2 q z ' 



(15) 



V(q, q z ,N c ) 



mq z 



2 2 

muj q 



muj 2 z q 2 z 



aN r h z 



Amqi 



2\l2-Kmq 2 q z 



nn. 



(16) 



We emphasise here that the form of the damping given by Eq. 1131 depends only on the relative angular frequency of 
the vortex and thermal cloud. This is a simple consequence of the presence of a thermodynamic instability of the 
vortex when < lj p . It should also be noted that the form that occurs in Ref. 01 is obtained from Eq. El by setting 
f2 = u) = a = 0, that is, only the kinetic energy contribution to cu p was included. 

Thus far we have retained the time dependence of the condensate occupation and phase, and it is pleasing to see 
that the vortex contributes a term to /z c (i) which depends on the relative rotation of the vortex and the cloud. To a 
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first approximation, the change in atom number during the vortex motion may be neglected, so that in what follows 
we fix the phase and occupation of the condensate. The time independent variational parameters q and q z are then 
chosen to minimise Eq. ED f° r N c atoms in the condensate. We choose to define the condensate length parameter q 
as the edge of the condensate for the vortex motion. 

For a non-rotating cloud the chemical potential decreases as the square of the distance of the vortex from the 
center of the trap, which is consistent with the thermodynamic instability of the vortex. If the cloud rotates at the 
frequency of the vortex precession in the laboratory frame, Q = lu p , the chemical potential becomes independent of 
the position of the vortex, and the vortex motion depends entirely on thermal fluctuations. When uj p < Q the center 
of the condensate becomes a minimum of the chemical potential, increasing quadratically with the radial distance of 
the vortex from the center of the trap. This is the region of thermodynamic stability. 



III. VORTEX DYNAMICS 



The influence of the thermal atoms on the motion of the vortex is the subject of the remainder of this paper, and 
we will use well established techniques to analyse the observable finite temperature effects on the vortex dynamics 
and its lifetime. 



A. Equations of motion 

As shown in Ref. , using techniques developed in Ref. |lOj , the equation of motion for the vortex coordinates can 
be mapped to a pair of equivalent Langevin field equations. Writing z(t) = u x (t) + iu y (t) for the vortex position in 
the complex plane, this takes the form of a single equation of motion 



dz 

lit 



(ioj p + j)z + r)(t) (17) 



where the complex Gaussian noise is determined by the non-zero correlator 

(v*(sHt)) = 2*6(8 -t), (18) 
and the magnitude of the thermal noise is determined by 



4N C 

For a vortex initially situated at z(0) = zq, the solution 



(19) 



z(t) = z a e^+^ t + / ds r ] {s)e^ + ^ t - s) , (20) 



together with Eq. ^| can be used to calculate all observables associated with the vortex motion. 

We can also find an equivalent Fokker-Planck equation (FPE) for the Langevin field equation derived above. Solving 
for the vortex motion is then equivalent to finding the probability distribution for the vortex position in the lab frame, 
which, from Eq. El can be seen to evolve according to the FPE 

d : (-iu p - j)z + -^{iu p - 7 )z* + ^T^zz) P(z t z*,t). (21) 



dt \dz dz* dzdz* 

This equation is closely related to the FPE for the P-function describing a quantised light field in a damped optical 
resonator |25j |. Indeed, the FPE is equivalent to the evolution of a damped Boson mode with frequency u) p , described 
by density matrix p, according to the master equation 

p = —iu> p \a" a, p] + (a — "f)(2apa i — a} ap — pa) a) 

+a(2a' pa — aa) p — paa'), (22) 

where [a, at] = 1, and a/ creates a particle in the oscillator mode with frequency u> p . The system is analogous to a 
one dimensional harmonic oscillator with a peculiar reservoir coupling, provided we regard u x and u y as conjugate 
phase space variables, analogous to the ordinary position and momentum. We will show below that the form of the 



6 



coupling leads naturally to the introduction of an associated quasimodc which neatly encapsulates the concept of 
vortex stability. The solution of Eq. 1211 for a vortex initially located at z is given by 

P(z,z ,t\z ,z ,0) = Me2yt _ 1) e XV ^ a{e2jt _ 1} j (23) 

which provides a complete description of the vortex motion for both thermodynamically stable and unstable regimes. 
The mean and variance of the distribution determine all higher order moments, and are given by 

(z(t)) = z e ( ^ +7)t , (24) 

(z*(t)z(t))-(z*(t))(z(t)) = -(e^-l). (25) 

7 

When the thermal cloud rotates slower than the vortex precesses the dissipation causes the mean radius to increase 
with time, indicating an instability, whereas if it rotates faster than the critical frequency it is damped with time, 
indicating stability. In both regimes thermal fluctuations cause the distribution to spread, so that the actual physical 
stability must be determined by including both processes. The sign of 7 determines the thermodynamic stability of 
the vortex, since this is a consequence of the sign of the precession frequency in the frame of the thermal cloud. As 
we will see below, thermal fluctuations weaken the thermodynamic stability, so that the vortex always has a finite 
lifetime. This is a consequence of diffusive instability which means that the the vortex can always reach the edge of 
the condensate in a finite time. 



B. Vortex Quasimode 



It is useful here to consider the angular momentum and energy associated with the vortex. We first consider the 
angular momentum corresponding to the Gaussian ansatz of Eq. El which is 

f ( ul(t) +ul(t)\ 

(L z ) = -ih / dx cf>*{x)(xd y - yd,)0(x) = hN c ( 1 - 2 yK ' , (26) 



so that 



/ u 2 Jt) +ul(t)\ 
N y {t) = N c I 1 - xW g2 yK> \ (27) 

is understood to be the number of particles contributing to the state with a vortex at radius [u^(t) + Uy(i)] 1 / 2 . We 
will also find it useful to introduce the number of particles associated with the motion of the vortex core, defined as 

n v (t)=N c -N v (t). (28) 

When the vortex is stationary at the center of the trap all atoms are in the vortex state, and the occupation of the 
mode associated with the vortex position is zero. 

To obtain an expression for the energy associated with the position of the vortex core we define the vortex Hamil- 
tonian 

H v = n v htu v (29) 

where 

uj v = Q — uj p (30) 

is the effective frequency of the mode. The average vortex energy is E v = (H v ), or 

E v (t) = (n v (t))huj v . (31) 

We can regard (n v (t)) as the mean occupation of the quasimode, and we will show that this interpretation is consistent 
with the equipartion theorem. For convenience we have chosen the origin of the vortex energy scale as Ml = HcOp, 
corresponding to the rotation rate of the cloud above which the vortex becomes thermodynamically stable. Note that 
the absence of a zero point energy term in Eq. 1291 is a consequence of working in the high temperature regime. 
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C. Vortex Lifetime I: Diffusion Time 



The vortex diffusion process is most easily characterised by the width the position distribution, ft is convenient 
to express the analysis of this section in terms of the rms radius r = \J (|z| 2 ), and we will use the notation N v (f) = 
N c (l — r 2 /q 2 ) to denote the mean occupation of the vortex state corresponding to the rms-radius f, and similarly for 
E v (f) and n v (f). We emphasise here that by using yj {\z\ 2 ) rather than (|z|) in our description of the lifetime we will 
avoid the need to introduce the approximation (l/|z|) w l/(|z|), which was implicitly used in the calculation of the 
vortex lifetime in Ref . [Tjj > effectively underestimating the role of fluctuations [2(| . 

We can now find an estimate for the vortex lifetime. The equation of motion for the rms radius 



f(f) - xj\zv\ 2 e 2 ^ + -\)h, 



(32) 



takes the form 



We can now find the time taken for f(t) to change from fi(t\) to ^2(^2) as 



(33) 



Using Eas. 11511151 and 01 we obtain 



r(ri,r 2 ) 



In 



In 



27 



1 - pE^n) 



l-pE v (f 2 ) 



cr/7 



cr/7 



(34) 



(35) 



In the above expression we have introduced r c representing the diffusion time obtained by choosing n = 0, r^ = q, 



and Q = u)p. which takes the form 



T(0,q) 



nhN r 



k B T J 24raa 2 ' 



(36) 



This universal critical diffusion time arises when the cloud rotates at the critical frequency for thermodynamic stability 
of the vortex, and its form is independent of the size of the condensate and the vortex precession frequency, arising 
only from thermal fluctuations which cause the vortex position to diffuse in time. This is the irreducible timescale 
associated with the diffusion instability. We note that the form of Eq. [JS] dictates that f\ < f-i in order to give a well 
defined positive value. 

The special case of experimental interest is the vortex lifetime associated with the time for a vortex initially located 
at the center of the trap to reach f = q. In Fig. we present this value for N c = 500 atoms of 87 Rb. As we noted 
above, the Gaussian ansatz is appropriate for weak interactions, which, for the species we consider, requires small 
condensate occupation. The expression for r c is shown by the dashed line, which also marks the boundary of the 
thcrmodynamically stable region, while the surface plot depicts the variation of Eq. 1351 with the temperature and 
angular frequency of the thermal cloud. 

To summarise, the diffusion time that we have computed in this section is expressed in terms of the variance of 
the vortex distribution (similar to the result of Ref. [H|). This approach reveals a characteristic timescale associated 
with the diffusive instability, and gives a good first estimate of the vortex lifetime. In the next section we give a more 
rigorous description of the vortex lifetime using the theory of mean first exit times [27l | . 



D. Vortex Lifetime II: Mean Exit Time 

We want to know the mean exit time associated with the escape of the vortex from the region r = {v^ + u 2 ) 1 / 2 < r%, 
given its initial position r = r%. After going into a rotating frame at frequency lu p , we may rewrite Eq. 1171 as the 
Langevin equation 

| =7^ (37) 
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FIG. 2: Vortex diffusion time for N c = 500 atoms of 87 Rb in an oblate trap with frequencies (ui z ,UJ r ) = 2ir x (10, 100)Hz. The 
diffusion time is shown for the stable range of Q {—uj r < fi < ui r ), and the bold line (color online) shows the critical lifetime r c . 



where the real Gaussian noise is defined by the correlation 

(C(*)C(*)>= 



(38) 



The symmetry of the system means that the vortex diffusion is an effectively one dimensional process in the radial 
coordinate. The problem of the mean exit time from the region (0, r 2 ) is formulated in terms of the probability that 
the vortex is in the region (0, r 2 ) at time t 



G(r u r 2 ,t)= / dr'p(r',t |n,0). (39) 
Jo 

Let the time that the particle leaves (0, 7*2 ) be T. Then Prob(T > t) = G(r\, r-x, t), and the mean first passage time is 



T(r u r 2 ) 



t dG(n,r 2 ,t). 



(40) 



The solution of this problem is known, and we refer the reader to Ref. [27| for the details [23|. The mean first exit 
time for the vortex initially at radius r\ to exit from the region r < r 2 is 



T(n,r 2 ) 



2pH v (q) 



In 



Bv(ri)e 



«l(-?Hv(r-l)) 



H v {r 2 )e E ^-P H ^2)) 



(41) 
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FIG. 3: Vortex diffusion time for N c = 800 atoms of 23 Na. The other parameters are the same as those used in Fig. [5] 



where we have made use of the exponential integral |29| , defined as 

r°° e"* 

E 1 (x)=J alt—. (42) 

Note that the singularity plays no role since the difference of these integrals is always over an interval which excludes 
the origin |3Cj . The mean exit time has similar order of magnitude to the expression found in section IlII CI given by 
Eq. |2S1 In particular, when the cloud and vortex co-rotate we have the critical exit time from the condensate center 
to the edge 

T c ^T(0,q) = §, (43) 

li — 10 v z 

which is a further manifestation of the diffusive instability. 

In Figs. 01 [3 we show the mean exit times for sodium and rubidium for a range of temperatures and angular 
frequencies of the thermal cloud. As before, we have marked the critical exit time T c with a dashed line. We note 
here that as well as being shorter by roughly a factor of three than the diffusion time estimate, the mean exit time 
has a weaker dependence on the relative motion of cloud and vortex than the diffusion time. 



10 




CO 

E 



q; 600 

o 



T[nK] 



^ [ co ] 

r J 



FIG. 4: Vortex mean exit time for N c = 500 atoms of 87 Rb for the same parameters as Fig. [5] 
E. Equipartition, Correlations and Fluctuations 

Using Eq. 1351 one can find an expression for the time dr required for the distribution to evolve from f to f + dr. 
After a little manipulation wc obtain the equivalent differential form 

^ = ^(l-^ v (r)), (44) 

which is just the rate equation for the vortex energy. The possibility of parametrising the energy by the radial diffusion 
time coordinate r arises from the physical fact that the mean energy associated with the vortex only depends on f . 
The solution is 

E v (t 2 ) = ^(r^e-^W^-n)/^ + keT ^ _ g-^fic^-n)/^ _ ^ 

When uj v > 0, corresponding to the thcrmodynamically stable scenario, we obtain the steady state occupation of the 
vortex quasimode 

(46) 



n v (r) 



hoj v ' 



in accordance with the equipartition theorem. As expected, the steady state distribution P s (z,z*), found from Eq. [231 
has the form 

N r huj v ( N c hui v \z\ 2 ' 



P s (z,z*) 



■ exp 



k B T q 2 y 



(47) 
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FIG. 5: Vortex mean exit time for iV c = 800 atoms of 23 Na for the same parameters as Fig. [5] 

which is consistent with the equipartition result of Eq. l4lfl 

Since we have such a simple description, we should also consider the measurable correlation functions for the 
vortex position that can give information about the thermal statistics of the condensate. It is convenient to discuss 
these quantities in the regime of thermodynamic stability, where two-time correlation functions reach a steady state 
behaviour. For the first and second order stationary two time correlation functions we have 



5 «(r) = lim t _ 



(z*(t)z(t + r)) 



9 {2) {r) 



lim. 



y/(z*(t) Z (t))(z*(t + T)z(t + T))' 
(z*(t)z*(t + T)z(t + T) Z (t)) 



(Z*(t)z(t))(z*(t + T)z(t + T))- 

The correlation functions associated with the steady state vortex motion are 



(48) 
(49) 



(50) 
(51) 



The Fourier transform of Eg. I5UI gives a Lorentzian power spectrum associated with the vortex postion. From Eg. 1511 
we see that there is a 'thermal bunching' effect [25( for the variance, corresponding to the increased probability of 
vortex- vortex coincidence at short times. Eq. E]also tells us about the steady state spectrum of energy fluctuations 
of the vortex quasimode, since we may rewrite this as 



(g v (0)g v (r))-(g v (0))(g v (r)) 
(H v (0))(H v (r)) 



-2] 7 |t 



(52) 
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giving a broad Lorentzian spectrum for the energy fluctuations arising from collisions with the thermal cloud. The 
power spectrum potentially has an important role in setting the thermal fluctuations arising in the condensate due to 
the vortex, since the vortex effectively acts as a radiating oscillator with this energy spectrum. However, the process 
of energy transmission into the condensate from a jostling vortex is expected to depend on the details of the vortex 
core structure, which is beyond the scope of this paper. 

Finally, we remark that the lifetime of the vortex can be estimated by short time measurements of the mean and 
variance of the vortex location. The radial drift rate is determined by 7 in Eg. 1241 while from Eg. 1251 we see that the 
radial variance evolves as Var(r) w 2at for short times, exactly as for the motion of a Brownian particle suspended in 
a viscous fluid. Since these two quantities, along with the vortex precession frequency and the angular frequency of 
the cloud, completely determine the vortex evolution, it is a simple matter to estimate the lifetime using Eg. 1351 once 
these are known. 



IV. CONCLUSION 

We have shown that a perturbative description which includes the rotation of the thermal cloud can be developed 
using the stochastic variational method of Duine et al [H,E3. The vortex acts as a probe for thermal fluctuations, with 
the lifetime depending on both the thermodynamic and diffusive instabilities through the temperature and motion of 
the thermal cloud. The process of vortex exit from the condensate can be satisfactorily described using the theory of 
mean exit times, and the theory predicts a shorter lifetime for the vortex than the diffusion time method. Nevertheless, 
the diffusive effect of the cloud is much stronger than the damping effect; in particular, if one neglects the noise from 
the description the exit time predictions are of the order of tens of diffusion times. 

The success of the approach relies on the physical fact that the variation of the atom number with vortex radius 
can be neglected to a first approximation. Even if the corrections arising from changing atom number were included, 
the method is not expected to give good predictions of the detailed motion near the condensate edge since it only 
affords a quadratic approximation to the chemical potential. Fortunately the motion in this region does not contribute 
significantly to the lifetime. Finall y, w e remark that the lifetime we have calculated bears little resemblance to the 
result of Fedichev and Shlyapnikov [31j , apart from a logarithmic behaviour which is characteristic of two dimensional 
systems. This is not unexpected since their analysis of vortex damping caused by quasiparticles in the Thomas- 
Fermi regime is for an entirely different scenario to the weakly interacting, thermally damped regime considered here. 
This raises the possibility of extending the high temperature classical field method beyond weak interactions to the 
Thomas-Fermi regime, which is a matter for future research. 
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